N91- 17721 I 


Theory of Plasma Contactors 
in Ground-Based Experiments and Low Earth Orbit 

M. J. Gerver* 

Plasma Fusion Center 
D. E. Hastings* and M. R. Oberhardt* 

Dept, of Aeronautics and Astronautics 
Massachusetts Institute of Technology 
Cambridge MA 02139 


Abstract 

Previous theoretical work on plasma contactors as current collectors has fallen into two 
categories: collisionless double layer theory (describing space charge limited contactor clouds) 
and collisional quasineutral theory. Ground based experiments at low current are well explained 
by double layer theory, but this theory does not scale well to power generation by electrodynamic 
tethers in space, since very high anode potentials are needed to draw a substantial ambient 
electron current across the magnetic field in the absence of collisions (or effective collisions 
due to turbulence). Isotropic quasineutral models of contactor clouds, extending over a region 
where the effective collision frequency v e exceeds the electron cyclotron frequency have 
low anode potentials, but would collect very little ambient electron current, much less than 
the emitted ion current. A new model is presented, for an anisotropic contactor cloud onented 
along the magnetic field, with v t < w«. The electron motion along the magnetic field is nearly 
collisionless, forming double layers in that direction, while across the magnetic field the electrons 
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diffuse collisionally and the potential profile is determined by quasineutrality. Using a simplified 
expression for i/ e due to ion acoustic turbulence, an analytic solution has been found for this 
model, which should be applicable to current collection in space. The anode potential is low 
and the collected ambient electron current can be several times the emitted ion current. 


1 Nomenclature 

Bo=ambient magnetic field 
c f =sound speed 

C= numerical factor relating electron thermal conductivity to electron transport 
e=charge on an electron 
£=electric field 

/ t =initial ionization fraction of source 
JV=azimuthal drag force on electrons 
<j=focussing factor due to anisotropy 
/ e =:e lectron current 
/,= ion current 
I = I t + /,=total current 

J e °°=ambient electron saturation current density 
J,*=ion current density 
fc x =perpendicular wave number 
L=Iength of tether 
m e =electron mass 

m, = ion mass 
n e —e lectron density 
n <a =ambient electron density 
n etf =cont actor electron density 

n, =ion density 

riic — contactor ion density 
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n, OU re*=80urce plasma density 
no=neutral density 
n 00 =electron density at infinity 
n cr i<.eai=neutral density required for ignition 
r==radial coordinate 
r a no<ie=anode radius 

r eore =radius at which electrons are collected, for any model 
»mn«r=i nner radius of double layer 
router— outer radius of double layer 
r lource =source radius 

ri=contactor cloud radius in anisotropic model 

Pload= load power 

i?,=tether resistance 

/J| oa j=load resistance 

T e£ =contactor electron temperature 

u*=electron flow velocity 

Uj=ion flow velocity 

u yt =Alfven speed 

uj=electron azimuthal drift velocity 
t; e = electron thermal velocity 
v r =radial velocity 

v # =axial velocity, parallel to the magnetic field 
vo=orbital velocity 
z=axial coordinate 

z 0 =half length of anisotropic contactor cloud 
/?«=ratio of electron pressure to magnetic pressure 
■floptical depth of source region to electron ionization 
A<£=potential drop across double layer 



q=electrical efficiency of the tether 
K=crosa field electron thermal conductivity 
Ap=Debye length 

Ax?,mn«r=Debye length at r inner 
Ao,ou<er=Debye length at r out<r 

/i=ratio of ion mass to proton mass ineffective electron collision frequency 

£ = ///<= gain 

p,=electron gyroradius 

<7=electron impact ionization cross-section 

^potential 

<^>o= anode potential 

<i>total — total tether potential 

c*> M =electron cyclotron frequency 

u pe =electron plasma frequency 


2 Introduction 

Plasma contactors are plasma clouds which allow the passage of charge between an electrode 
and an ambient plasma. They have been proposed for use in power generating devices such as 
electrodynamic tethers! *1 because they may substantially reduce the impedance of the electron 
current collection from the ionosphere and make the emission of electrons much less energetically 
expensive than using an electron gun. In this paper we will concentrate on plasma contactors 
used at an anode to collect electrons in the ionosphere or some other ambient plasma. Such a 
contactor will emit ions, as well as collect electrons. Two figures of merit for such a contactor are 
its impedance ^o//, and the gain £, defined as 

£ = 1 / A‘( r anod€)‘ 
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The impedance determines the maximum power that can be generated by a tether, since the total 
tether potential faotal ia fixed at v 0 BoL. If we ignore the ionospheric impedance and the impedance 
of electron emission, then 

$total ” Rioad^ -ft J H" 

The maximum power Ri oa d I 7 at fixed 4>totai and R% is obtained when R{ 0 ad = Rt + dfo/dl « 
Rt + <f>o/l. The power is greatest when the contactor impedance is lowest. The gain is important 
because it determines the rate at which gas must be used (to produce ions), for a given total current. 
If the gain is high, less gas is used to collect a given current. 

Both the impedance and the gain will depend on the current. In general there is a trade-off: at 
very low current, both high gain and low impedance are possible, but the power is low. While at 
high current, high gain can be obtained only at the cost of very high impedance (again resulting 
in low power). Low impedance and high power are possible only with low gain. To illustrate these 
trends, we may consider the extreme limits. When the current is equal to the electron saturation 
current of the ambient plasma over the surface area of the physical anode, then the gain is infinite 
(since no ions need be emitted to draw this much electron current) and the contactor impedence 
is zero, but the power (for low earth orbit and practical tether and anode parameters) is at most 
tens of watts. Arbitrarily large current (and high power) may be obtained by emitting a large 
ion current, but unless the anode potential is high enough, it will not be possible to collect many 
electrons across the magnetic field, and the gain will approach unity. A basic goal of contactor 
research is to determine how large a gain is possible at a given power level. If it turns out that 
at the power levels of interest for tethers (typically tens of kW) the maximum gain is close to 
unity, then there is no point in using plasma contactors for current collection; in effect, the best 
plasma contactor is no better than an ion beam. If, on the other hand, gains at least a few times 
greater than unity are possible at power levels of interest, then plasma contactors are useful as 
current collectors for tethers. We will present theoretical results suggesting that this is the case, 
although the gains are only moderate, in the range of 2 to 10. These theoretical results pertain 
to a regime (collisionless electron motion along the magnetic field, collisional diffusion across the 
magnetic field) which we expect to be valid in low earth orbit for high current contactors, but for 
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which there have been no ground based experiments. Such experiments are very important for 
confirming the theory, or showing how it must be modified. 

In previous work^^l it has been suggested that the plasma contactor cloud will consist of 
several different regions. First will be an inner core where the cloud will be isotropic because the 
two major directions of anistropy, namely the earth’s magnetic field and the direction of motion of 
the source will be shielded by the dense plasma from the contactor source. There will then be two 
outer regions where the two directions of anisotropy are manifested. Previously, it has generally 
been assumed that a substantial current of ambient electrons can be collected only from field lines 
that pass through the inner core regional. However, we will show in Section 4 that for conditions 
in low earth orbit it may also be possible to collect a significant electron current from the outer 
core region, where the anisotropy due to the magnetic field is important. 

There has been much debate about the size of the core region from which electrons can be 
collected. One estimate is obtained by matching the cloud density to the ambient density!®'®! 


n e( r core) ^ 

and another by taking magnetic field effects into account^) 

Vei'ccre) » Vce 

where v t is the radially dependent electron collision frequency (including effective “collisions” due 
to turbulence). A third estimate is obtained by requiring regularity of the self-consistent potential^ 

d *\ 

Or 'Tcor* 

and finally a fourth estimate comes by requiring a consistent space charge limited flow inside the 
core^l 

minni? \ Tcor , » m e n, uj|, €0r ,. 

These diverse theories give a wide range of current enhancement factors for the plasma cloud. 

If we assume a core cloud of radius r eore , then continuity of current gives 

f = Ii{ r <inod * ) + /« ( r anode) = 4( r c»») "b I*{ r cort) 
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and the gain is 


lejjcort ) 4(|tort) A( r ano<i<) j 
A { r arwd* ) A ( r anode) 

Plasma contactor clouds enhance or produce electron current flow through two possible paths. 
First (the first term on the right hand side of the equation), they can serve as virtual anodes 
through which electrons from far away can be drawn and collected to the real anode at the center 
of the cloud. Secondly (the second term on the right hand side), the neutral gas associated with 
the cloud can become ionized, creating electron-ion pairs. The electrons will be collected to the 
anode, and the ions will be repelled. For use in space with an electrodynamic tether, however, 
ionization of contactor neutrals is not an efficient use of neutral gas; if this is the only means by 
which the current is enhanced, then the same neutral gas can be used more efficiently by ionizing 
it internally in an ion source. Plasma contactors will be useful if they enable the ionosphere to 
supply electrons. The two sources of electrons in the ionosphere are the ionospheric plasma and the 
ionospheric neutrals. However the mean free path for ionization of the ionospheric neutral gas is so 
long (many kilometers) that ionization of this gas on the length scale of the plasma contactor cloud 
is highly unlikely. For this reason we shall assume that all ionization associated with contactors is 
ionization of contactor neutral gas. Therefore plasma contactors can be useful with electrodynamic 
tethers only if they enhance current by collecting ambient electrons from the ionosphere. The 
collected electron current A(r eore ) will generally be the saturation current times the area of the 
core cloud 4jrr^ t , or, if the contactor is only collecting electrons along magnetic field lines running 
into the core cloud, then 7 e (r eore ) will be the saturation current times 2xr 2 cort . (If, as we consider 
in Section 4, the core cloud is not spherical but is elongated in the direction of the magnetic field, 
then r eort is the minor radius, across the magnetic field.) For this reason the size of r eore is crucial 
to the effectiveness of plasma contactors as electron collectors in space. 

In Section 3 a collisionless double layer theory will be derived, along the lines of Wei and 
Wilburl 9 !, Amemiyal 10 l, and it will be shown that this theory provides a good quantitative de- 
scription of ground-based experiments at moderately low currents, but it will not be applicable to 
space-based contactors except at extremely low current and power. If the electrons are strictly col- 
lisionless, then the magnetic field prevents electrons from reaching the anode unless they originate 
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on field lines that pass close to the anode (which limits the current that can be collected) or the 
anode potential is high enough to pull electrons across the magnetic field to the anode from some 
distance away. A necessary condition for this, which depends on the anode radius r ano< ^, was found 
by Parker and Murphy^ Another constraint on r an0< k is that it must be less than the inner radius 
of the double layer. We will show that any spherically symmetric double layer with space-charge 
limited current greater than a very low limit (about 50 mA collected electron current, correspond- 
ing to 1 mA emitted ion current, for dayside equatorial low earth orbit, and even lower current for 
nightside) which satisfies these constraints must have an anode radius that is close to r cort . Such a 
plasma contactor would serve no purpose, since it would hardly collect any more ambient electron 
current than the bare anode. This means that an unmagnetized collisionless space-charge limited 
double layer model, as analyzed by Wei and Wilburt 9 !, cannot apply in space, except at very low 
currents, no matter how great the potential is. If the anode emits a current greater than this, at 
zero initial velocity (i.e. space-charge limited), and if the electrons are assumed to be collisionless, 
then the double layer cannot be spherically symmetric, regardless of the potential. Electron collec- 
tion will be inhibited across the magnetic field, and the collected electron current will be lower than 
predicted by the Wei-Wilbur theory^ for that anode potential and emitted ion current. Although 
a theory valid in this regime is not available, we can still obtain on upper limit on the collisionless 
electron current that can be collected, and a lower limit on the anode potential, for a given ion 
current, by assuming that the Parker- Murphy condition is marginally satisfied for a double layer 
obeying the equations of Wei and Wilbur, and ignoring the constraint that the inner radius of such 
a double layer must occur at a greater radius than r a no<u - We then obtain an upper limit to the 
power than can be generated by a plasma contactor collecting electrons to a 20 km long tether in 
space, in the absence of electron collisions. This maximum power is quite low, only a few hundred 
watts, less than an order of magnitude above the power that can be generated by a tether without 
a plasma contactor, using a bare anode to collect electrons. 

At higher emitted ion current, there will be a region where the electrons cannot go straight to 
the anode, but where ambient electrons will be trapped, to keep the plasma quasineutral. These 
electrons will rem ain trapped for a time long compared to the tune it would take for an unmagne- 
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tized electron to go straight to the anode. If there are effective collisions due to instabilities, some 
of these trapped electrons may be able to diffuse to the anode, and the collected electron current 
may be much greater than what would be found in the collisionless model. 

In Section 4, we will describe work on a model of the outer core region, in which the motion 
along the magnetic field is collisionless, forming a double layer, but the motion across the magnetic 
field is collisional and quasineutral. This model, which is expected to be applicable to contactors 
in space, suggests that significant current may be collected from this outer core region, with low 
contactor impedence. Unfortunately there are, to our knowledge, no experiments in this regime, 
to which the theory can be compared. Conclusions will be presented in Section 5 

3 Double-Layer Theory and Implications 

3.1 Collisionless Unmagnetized Model 

Ground-based experiments in which double layers are seen are well described by a collisionless 
unmagnetized model, as we will show. A schematic radial potential profile for such a model is 
shown in Fig. 1. We assume that the potential is monotonic, so there are two components of 
plasma, an ambient component and a contactor component. The ambient ions and electrons are 
maxwellian at positions r well beyond the double layer, with ion and electron temperatures T, a 
and T iat and density The contactor plasma has maxwellian electrons at temperature T tc and 
cold ions streaming radially out from a plasma source localized near the anode, with ion current 
/,*. The potential drop <j>o between the source, at r = r gourc4i and the ambient plasma at r — ► oo, is 
assumed to be much greater than any of the temperatures, and the radius at which the double layer 
forms is assumed to be much greater than a Debye length. With these assumptions, the plasma 
is quasineutral everywhere except inside the double layer, at r{ nntr < r < r 0U t €T . (Here r outer is 
the radius, called r cor4 in the Introduction, at which the ambient electron saturation current is 
collected.) Inside the contactor cloud, at r < r tnn€r , there are no ambient ions, and the density of 
ambient electrons, which have been accelerated in the double layer, is much less than the density of 
contactor electrons, so quasineutrality requires n <c (r) = n^r). The densities of contactor electrons 
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and ions are related to the potential <f> (defined relative to r — » oo) by 

n« e = n, ouree exp[(^ - )/T te ) (l) 

n, e = n ,ourct (’■<ou r e«/r) J [l + (<h - t)/T'c}~ l/2 (2) 

where we have assumed that ions are emerging from the source at the sound speed {T ec / m,) 1 ' ^ , 
due to acceleration in a Bohm presheath, and we have neglected any ionization or recombination 
occurring at r > r, ouree . Setting the right hand sides of Eqs. (1) and (2) equal to each other gives 
a transcendental equation for <j>(r). It is evident that for r » r (0ure( , 

^(r) » <£q — 2T ee ln(r/r source) (3) 

so the potential only drops a few times T tc inside the contactor cloud, much less than the total 
potential drop. The source density n J0urc< is related to the ion current /, by 

I* = 4xr* ource efi 90urC i (I^c/in,*) ^ (4) 


Outside the double layer, at r > r 0XiitTy the ambient electron density decreases from as r 
decreases, because no electrons are emitted from the double layer. We assume that there are no 
sources of electrons, or collisions, which could fill in the resulting empty region of velocity space. 
From quasineutrality, the ambient ion density must also decrease as r decreases (even if the density 
of contactor ions, accelerated in the double layer, is small compared to the ambient ion density), 
so the potential must rise by an amount on the order of Ti a . If T{ a is much less than T eay then 
the ambient electron density is not affected by the potential, so it is reduced from n <» by a simple 
geometric factor 

»..H = + (1 - (5) 


and the potential is given by 




( 6 ) 


(This rise in potential going from infinity to causes the ambient electrons to become supersonic 
by the time they reach fouttry so that they satisfy the Bohm sheath conditions’ll. This potential 
was calculated by Alpert, Gurevich and PitaevskiiSl for the case T{ a = T eay so we have labeled 
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this region the “Alpert-Gurevich presheath" in Fig. 1.) The potential drop from r outtr to oo is just 
T ia ln2, much less than the total potential drop. Most of the potential drop must therefore occur in 
the double layer. Within the double layer, r, nfM r < r < r outtr , the plasma is not quasineutral, and 
Poisson’s equation (for spherical symmetry) 

T'j/Tr = 4 ' ( "‘ ’ "•> < 7 > 

must be satisfied subject to the boundary conditions that <f> and d$fdr be continuous at Tinner and 
t outer* These four boundary conditions specify a solution to the second order differential equation, 
and the values of the free parameters r fflfier and r ou ttr * Since most of the drop in potential occurs 
in the double layer, to a good approximation the boundary conditions are 

&( r i nner) — ” 2T (ec ln(r,* fm€r /r JOUrc€ ) (8) 


pouter) — 0 


( 9 ) 


d<t>/dr - 0 at r inrur and r^er (10) 

If, as we have assumed, T, a T tat then the ambient ion density drops much more quickly than the 
ambient electron density as the potential starts to rise going inward from r outer , and we can neglect 
the ambient ion density in Eq. (7). Similarly, since the energy of the contactor ions is greater than 
T e e ri nn€r , even if only by a logarithmic factor, the contactor electron density drops much more 
quickly than the contactor ion density in going outward from r irm€r , and to rough approximation 
we can neglect the contactor electron density in the double layer. In the double layer, then, we 
must solve Poisson’s equation, Eq. (7), with 

n * = n ” - < ^ ,r exp(^/T ta )[l - erf(y^/T« a )] (11) 

'i o-Hr)Y 1/7 


n, = n. 


f r louree 


( 12 ) 


An approximate analytic solution, which provides some physical insight, may be found when the 
double layer is thin, i.e. »W <r - r, nneT < r inn(r . Then, in the vicinity of r inner , for A/j -C r - r inner <£ 


r outer — dinner, the potential approximates a Child- Langmuir sheath, with negligible n t 

\ 4/3 


4>( r inner) — ^(r) ft# ^ inner l r eource) | ^ 


inner 

.inner 


(13) 
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where 


( 14 ) 


is the Debye length at r,-. 


» j ^eelp( , *i’nr><r/^ SOi 

* D, inner ~~ 0-+2 « 

** c rijouree \ r source * 

, r . In the vicinity of r guter , for Xp ^ router — r r outer - r inner , the 


.(lin 

\r JO 


potential approximates an inverted Child- Langmuir sheath, with negligible n, 

*2rJ: -' 4/s 


* \ *D t outer ) 


where 


(15) 


(16) 


^D, outer — 2ir «l lla# 

is the Debye length at r^er- The transition from Eq. (13) to Eq. (15) occurs when n« ss n,-, which 
is to say at the point where the two expressions for 4>(r), Eqs. (13) and (15), have second derivatives 
that are equal in magnitude (but with opposite signs). At this point, the two expressions for ^(r) 
must have the same first derivative. This means that the transition from Eq. (13) to Eq. (15) 
must occur half way between r,' nn «r and r auter , with 4>(r) antisymmetric about this point, and the 
coefficients in front of the two expressions for must be equal, 

2T te ln(ri nntT /r,ource)X D J nntT — T ta XpJ )uier ( 1 7 ) 

Eq. (17) leads immediately to the well known double layer requirement! ? ! 


Ie/U = (18) 

where I e = 2xr* uter Jf 0 , and 7 e °° = en 0O (2xTt a /m e ) 1 ^ 1 is the ambient electron saturation current 
density. In other words, the contactor cloud will expand freely until the ion current density 7,/4xr* 
is equal to the ambient electron saturation current times (m t /rrii) 1 ^ . If T ta « T cc , then this will 
occur when the density of the contactor plasma is comparable to the density of the ambient plasma. 
From Eqs. (13), (15), and (17), the width of the double layer is related to the potential drop 

A (f> = ^{r inner) ~ 4>{ r outer) by 


r outer ~ r inner 


g outer 



(19) 


and these results are valid only if the width given by Eq. (19) is much less than r irmer . If this 
condition is not satisfied, then Poisson’s equation must be solved numerically, as has been done by 
Wei and Wilbur!®! and by Williams!*®!, and in this case /«//,• will be smaller than (: m,/m f ) 1 ^*. . 
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3.2 Comparison With Experiment 

The model outlined above is in good agreement with the ground-based experiments of Wilbur! 16 ' 

, in those conditions where double layers were seen. In these experiments, the anode had a radius 
r a nodt = 6 cm, but the effective source radius, where most of the ionization occurred, was r touree » 2 
cm. <fa could vary from 0 to 70V, and the collected electron current could vary from 0 to 1A. (At 
higher current, the effective collision frequency, due to streaming instabilities, was too high for 
collisionless double layer theory to be valid.) Neutral gas, xenon, was introduced at the center of 
the anode at a rate that could vary from 1.8 to 13.7 seem, which corresponded to a neutral density 
ranging from 3 x 10 11 to 10 1J cm' 3 , concentrated within »w« of the origin. For above some 
critical value, which depended on the neutral density, ambient electrons accelerated in the double 
layer had enough energy to ionize the gas, and the contactor cloud underwent a transition to an 
“ignited mode” in which this ionization was the major source of emitted ion current. The electron 
temperature and density and the plasma potential were measured as functions of position. The 
ambient ion temperature was much lower than the electron temperatures. 

In a typical case, with <h = 37V, most of the potential drop, 25V, occurred in a double layer 
(more or less spherical) located between r inntr = 8cm and r outer = 11cm. The rest of the potential 
drop occurred between the anode and r inner . The potential profile was virtually flat outside r outer . 
The ambient electron temperature was 5.5eV, and the ambient electron density was 3 x 10 7 cm 
These electrons have a Larmor radius of about 15 cm in the earth’s magnetic field, which is greater 
than r 0U t tT - r inntr , and once they cross the double layer they have a Larmor radius of about 50 
cm, which is greater than r* ut<r /2r an<><1 «, so the electrons can easily reach the anode according to 
the Parker-Murphy criterion! 11 !, and the assumption in our model of unmagnetized electrons is 
valid. The assumption of coUisionless electrons was also marginally satisfied if we estimate the 
effective collision frequency to be v t « 10- J w p «. At r ouler we find v t = 3 x 10V 1 , and the electron 
mean free path is about 30 cm, greater than the width of the double layer, while at r inntT we find 
j/ t = 2 x lO^s -1 and the mean free path of the accelerated ambient electrons is about 1 m, greater 
than Tinner- If the effective collision frequency is less than that taken here then the assumption of 
collisionless electrons is easily satisfied. Note that at densities a few times higher, the electron mean 
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free path would be comparable to the double layer width, and double layers could not exist. This 
is in agreement with observations at currents above 1A. There was also a 40eV ambient electron 
component (the “primary” electrons) of density 3 x 10 6 cm -5 . Such a component of electrons was not 
included in our model, but their effect can be included by using an effective T ea » 9eV which would 
give the same electron saturation current as that obtained from the 5.5eV and 40eV components. 

The collected electron current, 370mA, was in good agreement with this electron saturation 
current integrated over the area of the double layer 2*r£ jt4r (not 4 xr^ uter , since it was a half 
sphere). The electrons in the contactor cloud had a temperature T te = 2eV , and a density which 
went from 8 x 10 8 cm _s at r JOUr „ down to 2 x 10 7 cm -3 at r,-n», r . This ratio of n < (r toure( )/n ( (r tnn<r ) is 
close to {r inrur /r, <mret ) i [(<h - ^(r,nfi«r))/T„] 1/J , the value given by Eq.(2). The emitted ion current 
Ii would then be 2xr] ourct cn e {r, ouret ){T tc /m i y/ 2 = 0.4mA, fairly close to the ion current required 
by Eq.(17), (me/m,) 1 / 2 /, = 0.7mA. The observed width of the double layer, r outer - r inn «r « 3cm, is 
a few times greater than the width of 0.6 cm predicted by Eq.(18), but it is likely that the measured 
width is smeared out by fluctuations in the position of the double layer. Such fluctuations could 
be due to some intrinsic property of the double layer that would cause it to oscillate around 
equilibrium! instead of asymptotically approaching equilibrium. Such behavior is likely to be 
associated with non-monotonic potentials^®!, a feature that we have not included in our model. 
The fluctuations could also be caused by a more mundane effect, such as fluctuations in the gas 
feed. It would be of interest to try to measure such fluctuations and to determine their cause. 

3.3 L imi tations of Wei and Wilbur Model Due to Magnetized Electrons 

In Wilbur’s ground based experiments^®! the Larmor radius of the ambient electrons in the 
earth’s 0.3G magnetic field is about 20cm, much greater than the 3 cm thickness of the double 
layer, so the magnetic field will not significantly deflect the electrons as they cross the double layer. 
Once they cross the double layer, they will have a Larmor radius of about 50 cm, and in the 8 cm 
they have to traverse to get to the anode, they will be deflected by about j(8) 2 /50 = 0.7 cm, less 
than the 6 cm radius of the anode, consequently the magnetic field will not inhibit the electrons 
from getting to the anode!^!. Hence our model, which assumed unmagnetized electrons, ought 
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to be valid. An additional requirement of our model, r inntr > r anodt , is also satisfied in Wilbur’s 
experiments. 

In space, on the other hand, the ambient electron temperature, at least in the equatorial region, 
is much less, only about O.leV, so the Larmor radius is about 2.5cm, and the density is much less 
than in the ground based experiments (about 10 5 cm 3 rather than 3 x 10 7 cm 3 ). Therefore, to 
collect an electron current of several amps from the ambient plasma will require r outer of tens of 
meters, much greater than the electron Larmor radius. The electrons can traverse such a distance 
only if they undergo collisions (or effective collisions due to some kind of instability), or if they 
can gain enough energy as they cross the double layer to remain, in effect, unmagnetized, until 
they reach the anode. We have considered the latter possibility, and have found that, even with 
rather optimistic assumptions, it requires a sheath impedance that is undesirably large, since it 
would result in most of the tether potential drop occurring in the sheath. We conclude that effective 
collisions of some kind are needed in a plasma contactor in space, in order to collect a large electron 
current from the ambient plasma, at a reasonable impedance. 

Parker and Murphy! 11 ! have shown that, in the absence of collisions, and for e<fo » T eo , a 
necessary condition which must be satisfied for electrons at router to reach the anode is 

r luter/ r lnode < 1 + ( 8 tf< Ao/ m * w « r onod«) 1 ^ J ( 20 ) 

Eq. (20) is also a sufficient condition if all of the potential drop occurs in a thin double layer at r ouUr . 
If the double layer is thick, or if a significant part of the potential drop occurs in the quasineutral 
regions on either side of the double layer, then an even more stringent condition must be satisfied, in 
order for electrons to reach the anode. Another condition that must be satisfied is Tinner — T anodt • h 
turns out that for most parameters of interest this condition and Eq. (20) cannot both be satisfied, 
for a spherically symmetric space-charge limited collisionless double layer, as described by Wei and 
Wilbur! 9 ! and Amemiya! 19 !. This is true except at very low currents, or for anodes with r anodt 
almost equal to r outer- If higher ion currents are emitted from an anode (with rano<fe ^ r outer) with 
zero initial velocity, and there are no collisions or turbulence allowing electrons to be transported 
across the magnetic field, then a spherically symmetric double layer cannot develop, no matter how 
great the potential is. Electron collection will necessarily be inhibited in the direction across the 
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magnetic field; in this direction the potential profile will not follow the form found by Wei and 
Wilburl®!, because the collected ambient electron current will not be space-charge limited, but will 
be limited by magnetic field effects. A theory giving the electron current and potential in this 
anisotropic collisionless regime regime is not available. However, if we ignore the requirement that 
finntr > r ano<u and assume that only Eq. (20) and the Wei- Wilbur equations must be satisfied, then 
we can obtain an upper limit for the electron current than can be collected, and a lower limit for 
the potential, for a given ion current and anode radius. 

The electron current J« is related to r out „ by 

I* = ( 21 ) 

where J e °° = en ea {T ta /2‘xm t ) i ^ i is the ambient electron saturation current. We have calculated 
what the impedance of the double layer will be assuming Eq. (20) is barely satisfied, for r ano dt = 10 
cm. If, as turns out to be true, the resulting impedance is too high to make an efficient plasma 
contactor, we will know that we should look at plasma contactors in which the electrons undergo 
collisions (or are subject to turbulence which causes effective collisions) and diffuse into the anode, 
rather than going into the anode directly. 

Using Eq. (21) for I e , assuming Eq. (20) is barely satisfied, and using Wei and Wilbur’s 
calculation^) which relates r»ut«r Amrur uniquely to /*//,-, we can find and I t for a given 7, 
and electron saturation current J e °°. Since J e °° depends only on the properties of the ionosphere in 
low earth orbit, both I t and <£o are determined by /<. These values really represent an upper limit for 
I, and lower limit for 4>q, since Eq. (20) is only a necessary condition, not a sufficient condition, for 
collisionless electrons to reach the anode, and since we ignored the requirement that r inner > r anodt . 
The gain and potential drop are obtained by imposing the Parker-Murphy requirement and the 
limited source requirement (Eq. (21)) on the Wei and Wilbur results. 

In Figure 2 we show the gain ( against the ion current for argon and for a range of electron 
saturation current densities which span the range experienced in an equatorial low earth orbit 
(LEO). The gain is somewhat less than ( m,/m* )*/* = 272 for argon, and is weakly dependent on 
the ion current. We also show the associated potential drop through the double layer, which is 
really a lower limit on the potential drop. Typical potential drops are in the range of thousands 
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of volts for ion currents in the milliampere range. In Figures 2 through 4, the curves are dashed 
in the regime where Eq. (20) cannot be satisfied for a collisionless double layer with space charge 
limited current except by violating dinner ^ r anod * ■ Note that the curves are dashed except at the 


smallest ion currents, showing that a collisionless unmagnetized double layer with space charge 
limited current is not possible for most parameters of interest in low earth orbit. This conclusion 
does not depend on r an odf Making r anc j e < 10 cm would only make things worse, since, for a fixed 
ion current, r,- nn<r would shrink faster than r ano d*. Making ranod* much greater than 10 cm would 
allow higher ion and electron currents while satisfying Eq. (20) and > ranod * ■ However, for 

7 t °° < 2 x 10 -J A/m J , this could only be done if r an odt were nearly equal to r outer , in which case 
the plasma contactor would serve no purpose. Another way to show that this conclusion does not 
depend on ranod* i* to use Eqs. (19) and (20), with 4>o « A^, r ano d* — r,- nn<r , and r ou t er ss 2r irm<r . 
Combining these equations gives us 


*,-&(*)* 

(22) 

router « P * 

(23) 


where u> p * is the ambient electron plasma frequency and p* is the ambient electron gyroradius. 
Equations (18) and (21) then give 




(24) 


as the maximum ion current for which a coilisionless unmagnetized double layer with space-charge 
limited current is possible. This ion current depends only on ambient quantities and me/m,*, not 
on r ano< u or <£o, and is never greater than about 1mA for low earth orbit. 


In Figure 3 the total current is shown as a function of the electron saturation current density. 
The curve obtained for the collisionless double layer (really an upper limit) is shown for a fixed 
ion current of 10 mA. For comparison, we also show the total current for the isotropic quasineutral 
model described in Ref. [4], and for the anisotropic contactor model described in Section 4, for 
a fixed ion current of 1 Amp. This figure compares the realistic range of operation for the three 
models in typical ambient electron saturation current densities. A significant feature of this figure 
is that as the source varies by two orders of magnitude from 2 x 10“ 4 A/m J to 2 x 10~ 2 A/m 2 , the 
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total current (which is almost all collected electron current) varies by only a factor of 1.5, for the 
collisionless double layer model. This would seem to invalidate one of the conclusions in Ref. [l] 
which was that plasma contactors would not be useful on the nightside of an equatorial low earth 
orbit because the collected current would drop to almost nothing. Here the double layer moves 
out as the electron pressure drops so that the collected electron current is almost the same. On 
the other hand, if we took into account the actual requirements for electrons to reach the anode, 
rather than only using the Parker-Murphy condition, then it is likely that at low saturation current 
the double layer would be inhibited from moving out so far, and the collected electron current 
would be more sensitive to saturation current. Except for the upper end of the range of saturation 
current, the actual electron current that could be collected without collisions is certainly far less 
than the upper limit shown in Fig. 3. For the anisotropic collisional contactor model, which is 
more relevant for high current plasma contactors in low earth orbit, Fig. 3 shows that the total 
current is about 4 times higher, and the collected electron current is about 10 times higher, on the 
dayside (J e °° «2x 10“ 2 A/m 2 ) than on the nightside (Jf 0 «2x 1CT 4 A/m 2 ). 

In Figure 4 the current voltage characteristic is shown for the range of electron saturation 
current densities. At constant current in the milliampere range the voltage is seen to vary by two 
or three orders of magnitude for one order of magnitude variation in electron saturation current, for 
the collisionless double layer. At constant voltage, the current is roughly linear with the electron 
saturation current. Ampere range currents (which are mainly electrons) require tens of thousands 
of volts of potential drop, even for the highest value of the electron saturation current. These 
curves represent an upper limit on the electron current for a given potential, or a lower limit on 
the potential for a given electron current. For currents greater than about 50 mA, the space charge 
limited collisionless double layer model on which these curves are based cannot satisfy both Eq. (20) 
and ri nn€r > r amK ^; the actual potential needed to collect such currents, in the absence of collisions, 
would be far greater than the lower limits shown in Fig. 4. Curves for the isotropic quasineutral 
model and anisotropic model discussed in Ref. [4] and Section 4 are shown for comparison. 

With the use of these results we can calculate an upper limit for the current that could flow 
through a tether using a plasma contactor to collect electrons. A circuit diagram for a tether is 
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Table 1: Lo ad power against efficiency of double lay er contactor 


1 

Ii{m A) 

£ 

/(A) 

PloUW) 

0.1 

7 

26 

0.18 

100 

0.3 

6 

27 

0.16 

260 

0.5 

5 

35 

0.14 

380 

0.7 

2.5 

44 

0.11 

400 

0.9 

0.8 

75 

0.06 

290 


shown in Fig. 1 of Ref. [19]. The total potential drop fatal across the contactor, tether, load, 
and electron gun (or electron emitting contactor) is fixed by the length L of the tether, the earth’s 
magnetic field B 0 - 0.33 x 10 -4 T, and the orbital velocity of the spacecraft v 0 = 8km/s. For 
]j = 20km, we find <f> total = VoBqL — 5333V. The potential across the load is <j>ioad = Rioadi^i "1" !*)• 
The potential across the tether is Rt{I% + /«), where we take the tether impedance Rt = 2000. We 
could include the radiation impedance^®! but this is typically only about 10O, so may be neglected 
compared to the tether impedance. We also neglect the impedance of the electron gun or electron 
emitting contactor. If we assume a typical dayside ionosphere with J%° — 2 x 10 _J A/m : , a good 
fit to the numerical results in figure 4 is = &(/,• + 1*) 2 ' 08 where b — 1.8 x 10 s . For a given load 
Rioad, the current /=/,-- 1- 1, may be found by solving 

fatal = RtoadI+RtI+bI 20 * (25) 

and we may then find the power across the load Pioad — Rioad , and the efficiency r\ = Rtoad^I 4>totai , 
as functions of R loa d ■ (This definition of efficiency neglects the energy needed to produce the ions, 
but that is justified since this energy, about 50eV, is much less than the potential drop across 
the double layer, unless fj cs 99%. ) Table 1 shows Pioad and £ as functions of the efficiency 

n = Rioadt/fatar- 

The maximum power to the load is 400 W, but this occurs when the efficiency is only 70%. 
As noted in Ref. [l], in order for tethers to be competitive with other power systems in space it 
is necessary for them to operate at high efficiency, at least 80% or 90%. This is because all of 
the power has to be made up by periodically boosting the tether but only the load power can be 
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usefully used. If we desire an efficiency of 85%, then the maximum load power we can obtain is only 
320 W. The maximum power will in fact be much less than this, since Eq. (20) is not a sufficient 
condition for electrons to get across the magnetic field to the anode!* *1, and is known to be far 
from sufficient in the regime where r ou t tr r tnndr , which is true at the maximum power. Also, the 
requirement that r inner > r aru>< k is far from satisfied at the ion current needed for maximum power. 

We conclude that it is not possible to design a high power contactor which draws electrons 
straight across a double layer without collisions. Instead we should consider designs where collisions 
(or, more realistically, effective collisions due to instabilities of some kind) transport electrons across 
the magnetic field to the anode. 


3,4 Conditions for Ignited Plasma 

The calculations so far with the double layer model have all been for a totally ionized plasma. 
For a partially ionized plasma it is possible to include the effect of ionization and to show when the 
plasma will ignite. If we assume the neutral gas is expanding radially from the source at r toureti 
and that only a small fraction of it gets ionized, then the neutral density varies with radius as 
no(r) = no(r tource )(r iourct /r ) 2 . We apply conservation of mass from r tource to r, nrwr to obtain 


/eM = /e(w)exp(7(A^)[^ 


r »ourcc 

r inrur 


1 ) 


(26) 


where = no(r touret )r, ouret a. Here the electron ionization cross-section a is to be evaluated 

at a typical energy for an incoming ambient electron, &<f> + T ta . From conservation of current we 


obtain the gain as 


£=1 + 


- l)exp( 7 (l - ^“)) 


(27) 


where £(r, nrW r) = 
the double layer is 


1 + (£(»w) - 1)(1 - exp(7(l - ^“^))) 

). The ion current at the source in terms of the ion current just inside 


^4 = 1 + («w) - 1)(1 - «*PM1 - 7^))) (28) 

MA r mn«rJ “inner 

In order to interpret the calculations in Fig. 2 with ionization present we must interpret the ion 
current in the abscissa as /,(r irm€r ), and the gain as £(r infWf .). The relationship in terms of the ion 
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current emitted at the source is given above. It is apparent that there may be no positive solution 
of the source ion current for a given ion current at the double layer. Physically this will occur 
when there is so much neutral gas that the mixed gas-plasma flow ignites giving an avalanche of 
ion current. The ion current and collected electron current will continue to increase, and cannot 
reach a steady state until the collisionless double layer model is no longer valid. By setting the 
source ion current to zero we can obtain this critical neutral density for ignition as 

_ ~ ln(l ~ 1/ (( f inw<r)) 1 

^critical — j \ \^) 

V - 1 'source; r inn*r) T source** 

If we relate the source neutral density to the ion flow rate and initial fractional ionization (/,) we 
obtain ignition for 

fi ( r imwf) > Critic* (30) 


1 “A 


Taking r source = 0.1 m, c, = 4.89 x 10 3 m/s, a = = 3.21 x 1CT 20 m" 2 (for ionization 

of argon) and /* = 10~ 4 which is typical of hollow cathode devices, we find that the critical ion 
current is much greater, for a given gain, than the ion currents for which the collisionless double 
layer model is valid in low earth orbit, shown as solid curves occur in Fig. 2. Hence ignition will 
never occur in this regime. Ignition might be possible in the regime of higher ion current and lower 
gain typical of the anisotropic collisional contactor model described in Section 4. 


4 Anisotropic Contactor Model 

Hastings and Blandino^ considered a model where electrons were transported across the mag- 
netic field by effective collisions due to instabilities, and assumed that such transport could occur 
only out to a distance r C ore where the effective collision frequency v € was greater than the electron 
cyclotron frequency u/ Cf . With that model, they found that the collected ambient electron current 
for typical parameters in low earth orbit was less than the emitted ion current. Here, we consider 
the possibility that electrons can be collected from a more distant region where v t < u> ce . In that 
region the contactor cloud will be anisotropic, extending further in a direction along the magnetic 
field than across the magnetic field. We therefore use cylindrical coordinates z and r, where r now 
refers only to the distance across the magnetic field, not to the total distance from the anode as it 
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did in previous sections. We assume that the plasma density in the cloud is still great enough to 
short out the electric field due to the orbital velocity, so the cloud will be cylindrically symmetric. 
(At still larger distances from the anode, the effects of the orbital motion induced electric field 
will become important, and the cylindrical symmetry will be broken.) In this region the electron 
velocity will be mostly azimuthal, at the drift velocity 

= (M) 
m,uj et dr m e u ee dr m t w ee n t dr ' 

For parameters of interest, this drift velocity is much greater than the radial flow velocity of the 

emitted ions, which are effectively unmagnetized since we assume that the scale lengths are all 

much less than an ion Larmor radius. The velocity difference between the electrons and ions 

will then be nearly in the azimuthal direction. This relative cross-field drift velocity of magnetized 

electrons and unmagnetized ions can give rise to a several instabilities, among them the ion acoustic 

instability (both k±p t > 1 and k±p t < 1 varieties), the Buneman instability, the electron cyclotron 

drift instability (also known as the beam cyclotron instability), the modified two-stream instability, 

and the lower hybrid drift instability. Which of these instabilities dominates depends on such 

parameters as T«/T<, v^/c,, vj/v , , fj lt ui pt /ui ce , and vj/t/x- These instabilities will give rise to 

turbulent azimuthal electric fields, which will exert an azimuthal drag force Ft = v,mtv,i on the 

electrons, giving rise to a drift in the F x Bo (inward radial) direction at velocity 

Vr = —Vj ( 32 ) 

<*>«« 

We will assume that the potential drop in the plasma cloud is very much greater than the ion 
temperature Ti, which is typically only a few eV. Since, as we will show later, T t tends to be only a 
few times less than 4>q, this implies that T t /T{ 3> 1, except perhaps near the edge of the cloud. Also 
c$ •< V 4 ■< v t . In these circumstances, we expect the k ± p, > 1 ion acoustic instability to dominate 
(this is the same as the ion acoustic instability in an unmagnetized plasma). The effective collision 
frequency v, for this instability in its nonlinear saturated state scales with density like u pt , and is 
independent of for c, « t > 4 , but there is some uncertainty as to its dependence on T 4 /T< and 
v<t/ v »- We will simply assume that 

v t » lO -2 ^ (33) 
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independent of the other parameters. There is considerable theoretical and experimental evidence 
that the effective collision frequency due to ion acoustic turbulence is proportional to w p< > and 
somewhat weaker evidence that the constant of proportionality should be 10" J , as in Eq. (33). 
This evidence is discussed in Ref. [7). In addition, we note that particle simulations of saturated 
ion acoustic turbulence in infinite medium! 21 ’ 22 ' 23 ] generally give effective collision frequencies of 
this magnitude, and that experimental observations of collisionless shocks are in agreement with 
this result! 24 !, in a plasma contactor, the scale lengths are not infinite compared to the wavelengths 
of the unstable modes, the geometry differs from that of Ref. [24], and Eq. (33) may have to be 
modified. (Indeed, the requirement that the wavelengths of the dominant unstable modes be small 
compared to the radial scale length will probably set a lower limit on the ion current for which 
this model is valid.) A proper determination of v t would require a 3-D particle simulation of 
a contactor cloud, and experimental observations in the relevant regime to make sure that the 
simulation includes all of the relevant physics. Short of that, Eq. (33) is a reasonable guess that 
should be of help in choosing parameters for more careful theoretical and experimental studies. 
The method we will use to find analytic expressions for <j>(r,z ) and the collected electron current 
may also be applied using more realistic expressions for v t . 

The divergence of the radial flux of electrons due to v t and the radial electric field and temper- 
ature and density gradients must be balanced by an inward flux of electrons along the magnetic 


field, neglecting ionization and recombination: 


13 3 n 

m t v r + T~n t v t = 0 

r dr oz 


(34) 


At high densities, such as those in the experiment of Urrutia and Stenzel! 25 !, with a/ p , » u> ct , 
the mean free path of electrons will be short compared to the length of the contactor cloud, and the 
velocity v g along the magnetic field may also be found by balancing the force from the electric field 
ed<t>/dz with the drag force m,v,v,. In this case Eq. (34) will generally not be separable in r and 
z, and it is necessary to solve a fully two-dimensional partial differential equation. The boundary 
conditions will be that v r = 0 and <f> = 0 at the same surface, and the flux of electrons across this 
surface must be equal to the flux of the electron saturation current of the ambient plasma (along 
the magnetic field) outside the surface. The potential z) would be quasineutral everywhere. 
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Since the position of the ^ = 0 surface is not known in advance, this would be a difficult numerical 
problem. The ambient plasma in low earth orbit has much lower density, w,* < w ee , and this would 
also be true in most of a space-based contactor cloud, which, as we will show, would extend along 
the magnetic field to a distance where the cloud density is comparable to the ambient density. In 
this case, the electrons will flow freely along the magnetic field, and a different model is needed. If 
the total potential drop between the anode and the ambient plasma is greater than T te and T,a , 
then double layers will form at a distance z 0 along the magnetic field in both directions, where 

*-*($-&* (35) 
for thin double layers, just as in the unmagnetized collisionless case (see Eq. (18)). Here g(z) is a 
factor to take into account that the ions are focussed by the potential <£(r, z) if it is not spherically 
symmetric. Although the flow of electrons along the magnetic field is nearly collisionless, we will 
assume that there is enough drag to slow down the incoming electrons slightly, so that they will 
not escape out the other end, but will become trapped in the cloud. Only a small amount of 
drag is needed for this if cfo > T ea , and this could be provided by electron-electron streaming 
instabilities which produce effective collision frequencies of only a small fraction of w pe . Even if all 
of the electrons are not trapped, making this assumption will not introduce a large error if most of 
them are trapped. At z = ±zo, the flux of electrons along the field must then satisfy the boundary 
condition 

n e v, = (36) 

Because the flow of electrons across the magnetic field is collisional, no double layer exists in the 
radial direction. For fixed |z| < zq, <f>{r,z) must decrease smoothly to zero at some ri(z), satisfying 
quasineutrality all the way. For fixed r, along a given field line, as long as H r > z = 0) > T t (r), 
z) will not go to zero for |z| < zo. If is at least a few times greater than T„, then z = 0) 
will be greater than T t for all r not too close to ri(z = 0). It follows that rj is nearly independent 
of z. The contours of 4 > [ r t z )i snd the flow of ions and electrons, are shown schematically in Fig. 5. 
This means that Eq. (34) will be separable in r and z. The boundary conditions in r are 

Hr = r **04*, z) = <h + T e ln(n e (z)/n(z = 0)) (37) 
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at r = ri 


(38) 

(39) 


4>{r = ri) = 0 

d£ _ I dTe 

dr e dr 

The last condition follows from the fact that v r = 0 outside the contactor cloud, and there is no 
source or sink of electrons at r = rj, hence v r must vanish at ri just inside the contactor cloud. 
Eq. (31) (with T, = 0), and Eq. (32) then yield Eq. (39). 


4.1 Electron Temperature 


Before proceeding with the calculation of the potential profile <f>(r), we will briefly consider 
whether we are justified in assuming that is at least a few times greater than T t . The electron 
temperature profile T»{r) is determined by the balance between convection, conduction, and ohmic 
heating (both perpendicular and parallel to the magnetic field). We neglect ionization and line 
radiation, which should only be important near the anode, and we neglect heat lost by electrons 
boiling out along the magnetic field. 


-3 dT t 1 d dT t d<f> 

— v r — — H — rK— + ev r — 

2 dr rn, dr dr dr 


+ 



= 0 


(40) 


Here k is the cross-field thermal conductivity, which is dominated by turbulence just as the drag 
is. In general 


AC = 




(41) 


where C is a constant which depends on the details of the “collisions" causing the heat transport. 
For electron thermal conductivity across a magnetic field due to Coulomb collisions'^ , for example. 
C = 4.7. 

The boundary conditions are 


T, = 0 at r = ri (42) 

K a * — " ~ — n,v r T, at r = r ana< u (43) 

vf 4yr nyipjg zq 

where Q ia the heat flux going into the anode. This is generally greater than the convective heat 
flux into the anode (the second term on the right hand side), because (vj) for a half-maxwellian is 
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greater than (u*){v*). So dT t /dr > 0 at r ano j e . Because T t = 0 at r = rx, dT t /dr must change sign 
between r ano<it and rx, and we can estimate that the second term in Eq. (40) is of order -K.T t /n t r\. 
Using Eqs. (31), (32), and (41) we find 


v r = 


* ( dT e \ 

n t T t C \ dr dr ) 


(44) 


Then the first term in Eq. (40) is of order ±Ke<f>jCn t r\, and the third term is of order -t-«e 2 <£ 2 /Cn t T t r\ 
From Eqs. (34) and (36) the fourth term in Eq. (40) is comparable to (and has the same sign as) 
the third term. 

If C < 1 , it follows that the second and/or the first term must balance the third and fourth 
terms, so T t is of order e<f>. If C 3> 1, then the second term alone must balance the third and 
fourth terms, and T t » e^/C 1 / 2 «; e<f>. Our assumption that T t is at least a few times less than 4> 
is thus valid if C is somewhat greater than one. This is true for Coulomb collisions; whether it is 
true for ion acoustic turbulence is an open question that is beyond the scope of this paper. If k 
is dominated by an energetic tail of the electron distribution, perhaps electrons collected from the 
ambient plasma which have not yet thermalized, then C » 1. 


4.2 Potential Profile and Cloud Radius 

To find ^(r), we first integrate Eq. (34) over z from — zq to -t-zo, and use Eq. (36) to eliminate 

v, 

f+*0 1 d 

i ‘-rfr ’ n = 2J “ («> 

(It may seems counter intuitive that finding the electron current should require an integration over 
z, since no electrons are collected across field lines at the boundary of the cloud at r = rj, only 
along field lines at the boundary at z = ±zo. The purpose of the z integration is simply to show 
that the dominant contribution to the cross-field electron transport comes from small z, so that v t 
may be evaluated at z = 0 when the radial integration is done.) To obtain an expression for n«, 
which appears explicitly in Eq. (45) and also implicitly through the dependence of u t on u pt> we 
use quasineutrality 

", = n i = (4x)- 1 / < m | 1/1 e- 3 / 2 ( r 2 + z J )- l P (r,z)(*o - (46) 
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The expression for n,- in Eq. (46) comes from the fact that the ions are unmagnetized, and expanding 
spherically from the anode. The factor g{r, z) takes into account the focussing of the ions by 4>{r, z) 
which is not spherically symmetric. Using Eq. (32) for v r , Eq. (33) for u e , Eq. (46) for n ( , taking 
Bo = 0.3G, defining the ion atomic weight n = m,/m p , and expressing J, in amps, J e °° in amps/m l , 
and <f> and in volts, Eq. (45) becomes 


f 


+*> , 1 d 
dz - — 
r or 


r (^o - *r s/ v + ‘) d 4 = -i2/f s/ v s/ vr 


dr 


(47) 


Because (fo-^) and d<f>/dr are fairly independent of z, and the integrand is most strongly weighted 
near z = 0, we replace 4> and d<f>/dr by their values at z = 0, so they can be taken out of the integral. 
Similarly, we can set g{r,z) w 1, because self-consistently there cannot be a strong focussing effect 
for z < r where most of the contribution to the integral is. We then do the integration over z 


d_ 

dr 


> - *r v 'f r 


• = -12 r/,: 


-5/3 -3/4joo 


(48) 


We integrate Eq. (48) over r, using the boundary condition Eq. (39) to obtain the integration 
constant 

- *r sl, T- = e ir»V»/v.~(4 - >*) («) 


where 


,\ = r\+ irrv; s/ v>^i? / v /, W“)-‘ 


Finally we use Eq. (38) in Eq. (51) to obtain an equation for ri 

#J ' 4 = o .5 fn-WJT 


rt + ir 1 * 0 - ,/ v'^/? / y / vrr‘ 


(50) 


We integrate over r again, using Eq. (37) at z = 0 to obtain the integration constant 

_ ^)i/< = o.5/r 5/ V" 3/4 -/ e °°( 2r 2 rJ - f4 ) ( 51 ) 


(52) 


If, as we have been assuming, T, < efo, then the second term in brackets may be neglected, and 

ri = i .2^/“jfV /w (jrr 1/4 ( 53 ) 
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Note that r\ has an extremely weak dependence on <fc>. For almost any reasonable <£o, say 10V < 
o < 1000V, for argon, and for J e °° = 2 mA/m 2 , which is between the typical dayside and nightside 
values, 


r x « 15if /8 

(54) 

and 


/, = 2xrf « 2 lf f * 

(55) 

In general the total current /=/,■ + I € is 


/ = /, + 8(J~) 1/J 7? / V /, *o / * 

(56) 


A substantial ambient electron current can be collected for values of and total current that are of 
interest for tethers. For 1 A of argon at = 2mA/m 2 , for example, we get a gain ///, = 3, while 
for 0,5 A of xenon, at a typical dayside electron saturation current J € °° = 20m A/m 2 1 we obtain 
Ijli = 12. These gains, although not as large as the gains that were found with a completely 
collisionless double layer model, can still make a significant contribution to operation of tethers for 
power generation. These electron currents are much greater than the electron currents found in the 
quasineutral model of Ref. [4]; the physical reason for this is that electrons are transported across 
the magnetic field from much greater radius, where v t < u C€ . 

In Fig. 3, the total current is shown for a fixed ion current of 1 A, as a function of electron 
saturation current, using Eq. (56), and is compared to the total current for the isotropic quasineutral 
model discussed in Section 4, and for the collisionless double layer model using an ion current of 0.01 
A. Note that the current from Eq. (56) is much more sensitive to the electron saturation current 
than in the case of the collisionless double layer model. The reason is that the anisotropic contactor 
cloud, unlike the collisionless double layer cloud, cannot easily expand to larger radius to make up 
for a decrease in the ambient electron density. In Fig. 4, the current voltage characteristic is shown, 
from Eq. (56), for Jf° = 2 mA/m 2 , and compared to the results from the isotropic quasineutral 
model, and from the collisionless double layer model for a range of electron saturation currents. 
For realistic potentials, less than 1000V, the current from Eq. (56) is at least an order of magnitude 
greater than for the collisionless double layer model. 
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Table 2: Load power against efficiency of anisotropic contactor 



D 

m 

* 

/(A) 

PlndikW) 

0.1 

1 



2.64 

0.3 

1 

4.83 

4.83 

7.65 

0.5 

1 

4.6 

4.6 

12.1 

0.7 

1 

4.22 

4.22 

15.4 

0.9 

1 

2.32 

2.32 

10.6 


Table 3: Load power against efficiency of emitting an ion beam 



D 

m 

e 

/(A) 

Pload{kW) 

0.1 

22.64 

i 

22.64 

12.1 

0.3 

17.56 

i 

17.56 

28.1 

0.5 

12.48 

i 

12.48 

33.3 

0.7 

7.4 

i 

7.4 

27.7 

0.9 

2.3 

i 

2.3 

11.2 


Table 2 shows the load power F) 0 ad against efficiency, using the same ambient plasma and tether 
parameters as in Table 1, but using Eq. (56) to relate I and In this case, the maximum power 
obtained at ~ 80% efficiency is 12kW, much higher than in Table 1. Of course in a comparison 
with the collisionless double layer results the energetic cost of producing more ion current must be 
compared to the cost of the high potential associated with the space charge limited double layer. 

Finally in Table 3 we show the power to the load for a quasineutral model which just emits 
an ion beam, or a double layer model with ionization, so that a large current flows for very low 
potential drop (A^ ~ 0, £ = I). At 90% efficiency this configuration, which makes no use of 
the ambient plasma, can generate only slightly higher power than the anisotropic contactor, and 
requires substantially higher emitted ion current. This shows that the anisotropic contactor could 
make a significant contribution to the operation of tethers for power generation. 
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Table 4: Contactor models and the regimes where they are valid 


Model 

Limits of validity 

Applicable situations 

Where discussed 

Bare ion source 

/. > i. 

Ii > 1A in nightside LEO 

Table 3 

Bare anode 

u < u 

?outer ^ 2 r ano de 

Very low /<, large r anode , 
e.g. 1, < WfiA, r ano 4 t — lm, LEO 

Briefly in Sec. 2 

Collisionless 
unmagnetised 
double layer 

r outer > 2 f anode t r anode ^ r %nnert 
r outer / r anode < (c<^o/m e ) l ^UJ ce l , 
rLter ! r anode < No/ m e) * 

1/iA < U < 1mA in dayside LEO, 
ground-based experiments! 
with Ii < 1A 

Sec. 3 

Low 

collisionality 

magnetised 

r outer > 2 f anode i 
r outer / r anode > (*^o/ m «) 

u e < u f cet k ± r x < 1 

1mA < Ii < 100mA? in dayside LEO 

Not discussed 

Anisotropic 

(magnetised 

collisional) 

u e < u ee at r lt 
Pe < r u k ± r 1 > 1, 

r anode ^ **1 j 4 < Ie 

/, = 1A in dayside LEO 

Sec. 4 

Isotropic 

(unmagnetised 

collisional) 

v t > u ee at r u 
Vt/Ve < r lt 
> 1, Ii < I e 

Marginally in experiment 
of Urrutia tc Stensell^, 
never in LEO 

Briefly in Sec. 4, 
similar model 
in Ref.[4] 


5 Conclusions 

We have examined several models for electron collection by plasma contactors. The range of 
validity of the different models, and the situations where they are applicable, are summarized in 
Table 4. 

The ground based experiments at currents below 1 A appear to be well described by a double 
layer model which treats the electrons as collisionless and unmagnetized. In those experiments, 
the double layer forms approximately at the radius where the plasma emitted from the contactor 
reaches the ambient plasma density. This radius is less than or comparable to both the electron 
Larmor radius, and the mean free path of the electrons, based on a model for effective collisions due 
to instabilities. In high power space applications, where the plasma cloud must have a radius of 
tens of meters, and the ambient electron Larmor radius is only a few cm, neither of these conditions 
applies. Still neglecting collisions, but taking into account the finite electron Larmor radius, we 
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find that ambient electrons can get across the double layer and reach the anode only if the Parker- 
Murphy condition! 11 ! is satisfied (and even that is not a sufficient condition). For r ano<J< <5C r outer 
and ion current greater than the right hand side of Eq. (24) (about 1 mA for dayside low earth 
orbit, even lower for nightside), the Parker-Murphy condition cannot be satisfied for a spherically 
symmetric double layer with space charge limited current, since the r inntT determined by Wei and 
Wilbur! 9 ! would be less than r an ode, for any potential and r out „ satisfying the Parker-Murphy 
condition. This means that such collisionless double layers are not possible in space except at 
very low ion currents. This conclusion follows from the mass ratio and the magnetic field, electron 
density and temperature found in low earth orbit (since the right hand side of Eq. (24) depends only 
on these parameters), and does not depend on any assumption made about the potential or the size 
of the anode, other than r ano <u < r outer- (Collisionless double layers with higher ion currents are 
possible if r an0 ^ is made big enough so that the bare anode could collect almost as much electron 
current as the contactor cloud, but the contactor cloud would then serve no purpose.) At higher ion 
currents and small anodes, if we assume the electrons are still collisionless, the collected electron 
current will not be space charge limited, as assumed by Wei and Wilbur, but will be limited to a 
lower value by the magnetic field. Neglecting the requirement that r mn(r > r anoit , and considering 
only the Parker-Murphy condition, we found an upper limit to the collisionless electron current 
that could be collected, and a lower limit to the potential, as a function of ion current. We found 
that such a large potential is needed across the double layer in order to draw a reasonably large 
electron current that the available load power for a 20km long tether is never greater than 400 W . 
The maximum power is surely far less than this, since this figure was found for a configuration with 
*mfMr < r outer, and the Parker-Murphy condition is known to be far from sufficient in that limit; 
also, Tinner > r anode was known to be far from satisfied at the maximum power. The collisionless 
double layer model should be valid in space for emitted ion current sufficiently low (/, < 1 mA for 
dayside low earth orbit, much lower for nightside) that a double layer can form with 4>o < 5kV (the 
total tether voltage) allowing electrons to get across the magnetic field to the anode, and satisfying 
Unner > r anode- There is a further requirement for validity: the electrons must not be deflected 
from the anode by effective collisions, due to instabilities, as they are traversing the contactor. But 
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this requirement is easily satisfied in space, where the ambient u pe is not too much greater than 

w C€ . 

Since a plasma contactor described by the collisionless double layer model cannot generate 
anything close to the desired power, we must use much higher emitted ion currents. Although 
the transition from the collisionless double layer model to the collision&l quasineutral model is not 
completely understood, we expect at sufficiently high ion current that there will be instabilities 
strong enough to produce a high effective electron collision frequency in the contactor cloud. Such 
a contactor can be described by a colliaional quasineutral fluid model, in which electrons can flow 
across the magnetic field within a radius T core of the anode. If r corg is defined conservatively as 
the radius within which the effective electron collision frequency, due to ion acoustic and Buneman 
instabilities, exceeds the electron cyclotron frequency!^, then we find that the contactor has a 
very low impedance, but draws very little electron current because T cort is rather small. The total 
contactor current is hardly enhanced at all above the ion current that it is emitting. Even for 
those cases of higher T e where a modest gain in current occurs, that gain is due almost entirely to 
ionization of neutral gas emitted by the contactor, not to collection of electrons from the ambient 
plasma. In this case, the gas would probably be used more efficiently if it were ionized internally, 
in an ion source, rather than externally, where much of it can be lost. 

If we include the anisotropic part of the contactor cloud where the effective electron collision 
frequency is less than the electron cyclotron frequency, then electrons can be collected out to a 
much larger radius, and an electron current a few times greater than the ion current can be drawn 
from the ambient plasma, even at fairly low potentials. In contrast to the upper limits derived for 
the collisionless double layer model, and to the quasineutral model based on the more conservative 
definition of r cor «, the electron current has a significant dependence on the electron saturation 
current of the ambient plasma in this case, and is substantially higher, for a given ion current, 
on the dayside than on the nightside in equatorial low earth orbit. Analytic expressions for the 
potential profile and collected electron current can be obtained when the electron motion along 
the magnetic field is fairly coUisionless, so that a double layer forms in that direction, but the 
electrons flow collisionally across the magnetic field. This is the regime that is relevant to high 
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current plasma contactors in low earth orbit. Although the model which is solved analytically 
in Section 4 made the simple approximation that the effective electron collision frequency, due 
only to ion acoustic turbulence, is equal to 10 independent of T t and the electric field, 

the same method should be applicable using more realistic expressions for the effective collision 
frequency. Another approximation made in our analysis of this model is that there is sufficient 
electron thermal conductivity across the magnetic field to keep T t much lower than &> in the 
contactor cloud. The validity of this approximation must be examined using realistic turbulence 
models. If this approximation is at least marginally valid, then our results should be qualitatively 
correct. 

One important conclusion of our analysis is that most of the present ground based experiments 
have limited relevance to space applications of plasma contactors, since they operate in a regime 
where the magnetic field and effective collisions are not important, or only marginally important. 
This is true of space-based contactors only at very low current and power levels. An exception is 
the experiment of Urrutia and Stenzell 25 !, which examined a plasma in which the electron Larmor 
radius was small compared to the scale of the potential, and anomalous transport of electrons 
across the magnetic field was important. Indeed, they found that the anode collected an electron 
current a few times greater than the saturation current of the flux tube that intersected the anode, 
even when the effective collision frequency was less than the electron cyclotron frequency. Urrutia 
and Stenzel attributed their cross field electron transport to ion acoustic instabilities that were 
excited by the azimuthal E x B drift of the electrons relative to the unmagnetized ions, which 
gave rise to azimuthal wave electric fields which cause radial E x B drifts. In this respect the 
experiment was similar to the anisotropic contactor cloud model considered in Section 4. However, 
this experiment differed in one important respect from the regime, appropriate to low earth orbit, 
that was considered in Section 4. In the experiment, the density was about 2 x 10 xx cm -3 and 
w p «/w M « 50, much higher than in low earth orbit, and as a result the anomalous parallel resistivity, 
due to Buneman and ion acoustic instabilities excited by the relative electron and ion flow velocity 
along the field, was high. The electrons did not flow freely along the magnetic field, but diffused 
along the field like a collisional fluid, so there were no double layers along the field. It would be 
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desirable to do ground-based experiments in the regime where the electrons flow freely along the 
magnetic field but collisionally across the magnetic field, since this is applicable to high power 
plasma contactors in low earth orbit, and to compare the measured <f>(r y z) and collected current to 
the expressions calculated in Section 4, or to similar expressions found with more realistic models 
for v € . 

Another interesting feature seen by Urrutia and Stenzel is that the enhanced electron current 
was not continuous in time but occurred in periodic bursts, as the instabilities periodically grew up, 
saturated, and decayed. This behavior is probably due to the positive bias instability, which has 
been widely observed in configurations of this sort^I. It is not known whether similar behavior 
would occur in the regime of free electron flow along the magnetic field and collisional flow across the 
magnetic field, appropriate for low earth orbit. Theoretical and experimental studies are needed 
to answer this question, which could have important implications for power systems based on 
electrodynamic tethers in space. 


183 



References 


[1] M. Martinez-Sanchez and D. E. Hastings. A systems study of a 100 kW tether. J. of Astro- 
nautical Sciences, 35:75-96, 1987. 

[2] D. E. Hastings and N. A. Gatsonis. Plasma contactors for use with electrodynamic tethers for 
power generation. Acta Astronavtica, 17:827-836, 1988. 

[3] E. P. Szuszczewicz. Technical issues in the conduct of large space platform experiments in 
plasma physics and geoplasma sciences. In Space Technology Plasma Issues in 2001, Jet 
Propulsion Laboratory 86-49, California Institute of Technology, Pasadena, California, 1986. 

[4] D. E. Hastings and J. Blandino. Bounds on current collection by plasma clouds from the 
ionosphere. Journal of Geophysical Research, 94:2737-2744, 1989. 

[5] D. Parks, M. J. Mandell, and I. Katz. Fluid model of plasma outside a hollow cathode 
neutralizer. Journal of Spacecraft and Rockets, 19:354-357, 1982. 

[6] D. Parks and I. Katz. Theory of plasma contactors for electrodynamic tethered satellite 
systems. Journal of Spacecraft and Rockets, 24:245-249, 1987. 

[7] D. E. Hastings. Theory of plasma contactors used in the ionosphere. Journal of Spacecraft 
and Rockets, 24:250-256, 1987. 

[8] M. Dobrowolny and L. less. Model of the interaction of a hollow cathode with the ionosphere. 
In Second International Conference on Tethers in Space, American Institute of Aeronautics 
and Astronautics, 1988. 

[9] R. Wei and P. Wilbur. Space charge limited current flow in a spherical double sheath. Journal 
of Applied Physics, 60:2280-2284, 1986. 

(10] H. Amemiya. Space charge limited current between coaxial cylinders and concentric spheres. 
Scientific Papers of the Institute of Physical and Chemical Research, 63:1—6, 1969. 



[11] L. W. Parker and B. L. Murphy. Potential buildup of an electron emitting ionospheric satellite. 
J. Geophys. Rea., 72:1631, 1967. 

[12] Lewi Tonks and Irving Langmuir. A general theory of the plasma of an arc. Physical Review 
34:876, 1929. 

[13] David Bohm. Minimum ion kinetic energy for a stable sheath. In The Characteristics of 
Electrical Discharges in Magnetic Fields (A. Guthrie and R. K. Wakerling, eds.), McGraw- 
Hill, 1949. 

[14] Ya. L. Al’pert, A. V. Gurevich, and L. P. Pitaevskii. Space Physics with Artificial Satellites. 
Consultants Bureau, 1965. 

[15] J. D. Williams, private communication. 

[16] P. Wilbur and T. Laupa. Plasma contactor design for electrodynamic tether applications. 
Advances in Space Research, 8:221-224, 1988. 

[17] M. A. Raadu and M. B. Silevitch. On the negative resistance of double layers. J. Plasma 
Phys., 30:249-254, 1983, and references therein. 

[18] K. Y. Kim. Theory of nonmonotonic double layers. Phys. Fluids, 30:3686-3694, 1987. 

[19] J. D. Williams and P. J. Wilbur. Plasma contacting- an enabling technology. 27th Aerospace 
Science Meeting, January 9-12, 1989, Reno, Nevada, AIAA-89-0677. 

[20] D. E. Hastings and J. Wang. The radiation impedance of a electrodynamic tether with end 
connectors. Geophysical Research Letters, 14:519-522, 1987. 

[21] D. Biskamp and R. Chodura. Computer simulation of anomalous dc resistivity Phys. Rev. 
Lett., 27:1553-1556, 1971. 

[22] M. Z. Caponi and R. C. Davidson. Influence of ion-resonance broadening on the anomalous 
heating and momentum transfer in a current-carrying plasma. Phys. Rev. Lett., 31:86, 1973. 

[23] J. A. Wesson and A. Sykes. Theory of ion-sound resistivity. Phys. Rev. Lett., 31:449, 1973. 


185 



[24] P. C. Liewer and N. A. Krall. Self-consistent approach to anomalous resistivity applied to 
theta pinch experiment. Phys. Fluids, 16:1953-1963, 1973. 

[25] R. L. Stenzel and J. M. Urrutia. Labratory model of a tethered balloon-electron beam current 
system. Geophysical Research Letters, 13:797—800, 1986. J. M. Urrutia and R. L. Stenzel. 
Anomalous currents to an electrode in a magnetoplasma. Phys. Rev. Lett. 57:715, 1986. 

[26] David L. Book. NRL Plasma Formulary. Naval Research Laboratory, 1980. 




r«odi r inn*f r out«r 


Figure 1 Schematic radial potential profile for collisionless unmagnetized double layer. The 
Bohm presheath is described in Ref. [12] and [13], and the Alpert-Gurevich presheath 
in Ref. [14]. 



ION CURRENT IA.p.r.,1 

Figure 2 Gain and lower limit on potential drop, as functions of the emitted argon ion current 
and the electron saturation current, for collisionless double layer with space charge 
limited current, marginally satisfying the Parker-Murphy condition with a 10 cm anode 
radius, using Eq. (20) and Eq. (21), and Fig. 5 of Ref. [9]. 
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Figure 3 Total collected current vs. electron saturation current with the emitted ion cur- 
rent held constant, for the collisionless double layer model (upper limit), the isotropic 
quasineutral model^l and the anisotropic contactor model. 



Figure 4 Total current vs. potential drop for the collisionless double layer model, the isotropic 
quasineutral model^l and the anisotropic contactor model. 
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Anisotropic Contactor Model 



Figure 5 Schematic picture of the anisotropic contactor modei t showing equipotential contours 
and the flow of ions and electrons. 


